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Abstract 


The  grid  sensitivity  of  the  Reynolds-averaged  Navier-Stokes  (RANS)  solvers  ANSYS 
CFX  (ANSYS,  Inc.)  and  TRANSOM  (DRDC  Atlantic)  has  been  tested  for  two- 
dimensional  flows  around  propeller  blade  sections.  This  is  a  preliminary  step  before 
developing  grids  for  the  solution  of  flows  around  cavitating  propellers.  The  sensitivity 
of  the  pressure  and  skin  friction  distributions  near  the  leading  edge  were  tested  for 
variations  in  near  wall  node  spacing,  the  node  density  near  the  leading  edge,  and  the 
overall  size  of  the  grid. 


Resume 


La  sensibilite  du  maillage  des  solutionneurs  de  1’ Analyse  d’equations  de  Navier-Stokes 
moyennees  (RANS)  ANSYS  CFX  (ANSYS,  Inc.)  et  de  TRANSOM  (RDDC  Atlan- 
tique)  a  ete  eprouvee  pour  les  ecoulements  bi-dimensionnels  autour  de  sections  de 
pales  d’helice.  II  s’agit  d’une  etape  preliminaire  avant  I’elaboration  de  maillages  pour 
la  solution  des  ecoulements  autour  des  helices  avec  cavitation.  La  sensibilite  des  dis¬ 
tributions  de  pression  et  de  resistance  de  friction  pres  du  bord  d’attaque  ont  ete 
eprouvees  en  vue  de  deceler  les  variations  dans  I’espacement  de  noeud  pres  des  parois, 
la  densite  de  noeud  pres  du  bord  d’attaque  et  les  dimensions  hors-tout  du  maillage. 
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Executive  summary 

Grid  dependence  of  RANS  codes  for  flows  past 
propeller  blade  sections 

David  Hally;  DRDC  Atlantic  TM  2008-262;  Defence  R&D  Canada  -  Atlantic;  April 
2009. 

Background:  The  Maritime  Asset  Protection  Section  at  DRDC  Atlantic  is  applying 
Reynolds-averaged  Navier-Stokes  (RANS)  solvers  to  the  prediction  of  the  flow  around 
propellers  so  that  the  noise  generated  by  propeller  cavitation  can  be  predicted.  An 
important  part  of  this  research  effort  is  to  determine  how  best  to  generate  the  compu¬ 
tational  grids  around  the  propeller.  Since  propeller  blades  are  made  from  a  series  of 
airfoil-shaped  sections,  as  a  preliminary  step,  the  grid  sensitivity  of  two-dimensional 
RANS  calculations  around  propeller  blade  sections  was  addressed  first. 

Principal  results:  The  sensitivity  of  the  RANS  solvers  ANSYS  CFX  (ANSYS,  Inc.) 
and  TRANSOM  (DRDC  Atlantic)  to  variations  in  computational  grids  has  been 
assessed  for  flows  around  propeller  blade  sections.  The  sensitivity  of  the  pressure  and 
skin  friction  distributions  near  the  leading  edge  were  tested  for  variations  in  near  wall 
node  spacing,  the  node  density  near  the  leading  edge,  and  the  overall  size  of  the  grid. 

Significance  of  results:  These  calculations  provide  information  on  how  computa¬ 
tional  grids  on  propellers  should  be  constructed  so  that  accurate  predictions  of  the 
pressures  distribution  can  be  obtained.  This  is  a  preliminary  step  in  an  attempt  to 
calculate  the  flow  over  full  propellers  so  that  cavitation,  and  ultimately  radiated  noise 
levels,  can  be  predicted. 

Future  work:  Procedures  for  generating  grids  on  propellers  will  be  developed.  The 
grids  will  be  used  with  ANSYS  CFX  and  TRANSOM  to  calculate  flows  past  propellers 
for  which  good  experimental  data  are  available.  Once  accurate  predictions  of  the  flow 
are  attained,  procedures  for  calculating  the  propeller  cavitation  will  be  developed  and 
used  to  develop  a  radiated  noise  model. 
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Grid  dependence  of  RANS  codes  for  flows  past 
propeller  blade  sections 

David  Hally ;  DRDC  Atlantic  TM  2008-262 ;  R  &  D  pour  la  defense  Canada  - 

Atlantique ;  avril  2009. 

Contexte  :  La  section  de  Protection  des  biens  maritimes  de  RDDC  Atlantique  ap¬ 
plique  les  solutionneurs  de  1’ Analyse  d’equations  de  Navier-Stokes  moyennees  (RANS) 
a  la  prediction  de  Lecoulement  autour  des  helices  afin  de  pouvoir  predire  le  bruit  pro- 
duit  par  la  cavitation  de  I’helice.  Get  effort  de  recherche  vise,  en  grande  partie,  a 
determiner  comment  produire  les  meilleurs  maillages  informatiques  autour  de  I’he¬ 
lice.  Puisque  les  pales  d’helices  sont  fabriquees  d’une  serie  de  sections  en  forme  de 
proffl  aerodynamique,  en  tant  qu’etape  preliminaire,  la  sensibilite  de  maillage  des 
calculs  RANS  bi-dimensionnels  autour  des  sections  de  pales  d’helice  a  ete  abordee  en 
premier. 

Resultats  :  La  sensibilite  des  solutionneurs  RANS  ANSYS  CFX  (ANSYS,  Inc.)  et 
TRANSOM  (RDDC  Atlantique)  aux  variations  des  maillages  informatiques  a  ete 
evaluee  pour  les  ecoulements  autour  de  sections  de  pales  d’helice.  La  sensibilite  des 
distributions  de  pression  et  de  resistance  de  friction  pres  du  bord  d’attaque  a  ete 
eprouvee  en  vue  de  deceler  les  variations  dans  I’espacement  de  noeud  pres  des  parois, 
la  densite  de  noeud  pres  du  bord  d’attaque  et  les  dimensions  hors-tout  du  maillage. 

Importance  :  Ces  calculs  fournissent  de  I’information  sur  la  fagon  dont  les  maillages 
informatiques  sur  les  helices  doivent  etre  fabriquees  afin  de  pouvoir  obtenir  des  pre¬ 
dictions  precises  de  la  distribution  des  pressions.  II  s’agit  d’une  etape  preliminaire 
dans  une  tentative  de  calcul  de  I’ecoulement  sur  des  helices  completes  de  fagon  a 
pouvoir  calculer  les  niveaux  de  bruit  de  cavitation  et,  a  la  limite,  de  bruit  rayonne. 

Perspectives  :  Des  methodes  de  production  de  maillages  sur  helices  seront  elaborees. 
Les  maillages  seront  utilises  avec  I’ANSYS  CFX  et  TRANSOM  pour  calculer  les 
ecoulements  sur  les  helices  pour  lesquels  de  bonnes  donnees  experimentales  sont  deja 
disponibles.  Une  fois  que  des  predictions  precises  de  I’ecoulement  seront  obtenues,  des 
methodes  pour  le  calcul  de  la  cavitation  d’helice  seront  elaborees  et  utilisees  pour  le 
developpement  d’un  modele  de  bruit  rayonne. 
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1  Introduction 


The  Maritime  Asset  Protection  Section  at  DRDC  Atlantic  wishes  to  apply  Reynolds- 
averaged  Navier-Stokes  (RANS)  solvers  to  the  prediction  of  cavitation  on  marine 
propellers.  One  question  of  interest  is  how  dense  the  grid  must  be  on  the  blade 
surface  in  order  to  model  the  pressure  adequately.  The  answer  will  have  important 
implications  for  the  gridding  strategy  used.  Since  over  most  of  a  propeller  blade 
the  flow  can  be  approximated  by  two-dimensional  flow  over  a  series  of  airfoil-shaped 
sections,  the  simpler  question  of  the  required  grid  density  for  airfoils  was  tackled  first. 

The  DTMB  modification  [1]  to  the  NACA  66  airfoil  was  used  for  this  study.  It  has 
often  been  used  as  a  section  for  propellers. 

A  series  of  grids  were  generated  using  the  grid  generator  Pointwise  [2] .  Scripts  were 
written  using  the  Pointwise  scripting  language  Pointwise  Glyph  (based  on  the  pro¬ 
gramming  language  Tcl)  to  allow  the  automatic  generation  of  grids  suitable  for  use 
with  the  RANS  solvers  currently  used  by  DRDC  Atlantic:  ANSYS  CFX  [3],  a  commer¬ 
cial  code  from  ANSYS,  Inc.;  and  TRANSOM,  a  code  developed  at  DRDC  Atlantic  [4]. 
Each  grid  depends  on  a  number  of  parameters  which  were  then  varied  to  determine 
their  effects  on  the  pressure  and  skin-friction  distributions.  Particular  importance 
was  placed  on  the  flow  near  the  leading  edge  where  cavitation  will  normally  appear. 

It  is  important  to  note  that  this  study  was  neither  a  verihcation  nor  a  validation 
exercise.  There  are  no  comparisons  with  experiment,  nor  with  theoretical  solutions. 
Both  ANSYS  CFX  and  TRANSOM  have  been  validated  extensively  [5,6].  It  is  merely 
a  study  to  determine  what  grid  density  is  necessary  to  produce  results  that  do  not 
depend  on  the  grid. 

2  Procedure  for  making  grids  in  Pointwise 


Two  scripts  were  written  in  Pointwise  Glyph:  one  to  generate  grids  for  ANSYS  CFX 
and  the  other  to  generate  grids  for  TRANSOM.  Each  starts  with  the  geometry  of  the 
airfoil  in  ICES  [7]  format. 

An  0-topology  was  chosen  for  the  grids  as  this  is  the  topology  planned  for  gridding 
propellers.  The  trailing  edge  of  the  airfoil  is  closed  using  a  nearly  semi-circular  exten¬ 
sion:  see  Figure  1.  This  makes  it  easier  to  implement  the  0-topology.  On  a  propeller, 
a  similar  closure  at  the  trailing  edge  causes  the  whole  propeller  blade  to  be  smooth 
with  no  mismatch  at  the  tip  between  a  smooth  leading  edge  and  sharp  trailing  edge. 
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Figure  1:  A  trailing  edge  before  (left)  and  after  (right)  closure  with 
a  semi-circular  extension  (right). 


2.1  The  TRANSOM  grid 

The  TRANSOM  grid  is  structured  and  uses  0-topology.  It  is  generated  by  extrusion 
from  the  airfoil  surface.  The  grid  is  controlled  by  six  parameters  each  of  which  is 
given  a  nominal  value.  The  parameters,  with  their  nominal  values  in  brackets  (all 
distances  are  relative  to  the  chord  length),  are: 


1.  The  mean  spacing  of  the  nodes  on  the  airfoil  (0.02). 

2.  The  number  of  nodes  per  radius  of  curvature  at  the  leading  edge  =  5). 

3.  The  number  of  nodes  per  radius  of  curvature  at  the  trailing  edge  {Nte  =  5). 

4.  The  node  spacing  between  nodes  on  the  airfoil  and  those  just  off  it  (10“^). 

5.  The  rate  at  which  node  spacing  grows  close  to  the  airfoil  (1.1). 

6.  The  approximate  outer  radius  of  the  grid  (i?  =  20). 

The  grid  is  constructed  as  follows: 


1.  The  spacing  of  the  nodes  at  the  leading  edge  and  the  trailing  edge  is  determined 
according  to  the  values  of  A)e  and  Nte  and  the  thickness  of  the  airfoil.  (The 
leading  edge  radius  of  curvature  varies  as  the  square  of  the  thickness,  the  trailing 
edge  curvature  in  proportion  to  the  thickness.) 


2. 


The  node  spacing  at  the  leading  edge  is  allowed  to  increase  geometrically  by  a 
factor  of  1.1  nntil  the  spacing  between  nodes  is  equal  to  the  mean  spacing.  This 
reqnires 


N 


log{Sm/sie) 

log(l.l) 


+  2 


(1) 
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Figure  2:  An  example  of  a  grid  used  by  TRANSOM.  The  full  grid  is  shown  on  the 
left,  a  close-up  near  the  airfoil  on  the  right.  This  grid  is  coarser  near  the  leading  and 
trailing  edges  than  the  grids  actually  used. 


nodes  on  each  side  of  the  leading  edge,  where  Sm  is  the  mean  spacing,  s/e  is  the 
spacing  at  the  leading  edge  and  |a;|  denotes  the  closest  integer  to  x. 

3.  The  nodes  near  the  trailing  edge  are  handled  in  a  similar  fashion  to  those  at 
the  leading  edge. 

4.  The  remainder  of  the  airfoil  is  covered  by  nodes  whose  spacing  is  the  mean 
spacing. 

5.  The  nodes  at  the  airfoil  surface  are  extruded  for  steps  using  the  Pointwise 
hyperbolic  extrusion  with 


log(l  +  R{w  -  l)/g^) 
log(u;) 


(2) 


where  R  is  the  outer  radius,  w  is  the  extrusion  growth  rate  and  s^j  is  the  node 
separation  at  the  wall.  A  Kinsey-Barth  smoothing  factor  of  3  is  used  during 
the  extrusion. 


6.  The  grid  is  saved  in  PLOT3D  format. 


An  example  of  a  TRANSOM  grid  is  shown  in  Figure  2. 
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2.2  The  grid  for  ANSYS  CFX 

The  grids  used  with  ANSYS  CFX  are  hybrid  with  an  inner  structured  block  with 
0-topology  and  an  outer  unstructured  block.  The  outer  boundary  is  a  square  (this 
simplihes  the  setting  of  boundary  conditions  in  the  ANSYS  CFX  pre-processor).  The 
inner  block  is  generated  by  extrusion  from  the  airfoil  surface.  The  grid  is  controlled 
by  nine  parameters  each  of  which  is  given  a  nominal  value.  The  parameters,  with 
their  nominal  values  in  brackets  (all  distances  are  relative  to  the  chord  length),  are: 

1.  The  mean  spacing  of  the  nodes  on  the  airfoil  (0.02). 

2.  The  number  of  nodes  per  radius  of  curvature  at  the  leading  edge  {Nie  =  5). 

3.  The  number  of  nodes  per  radius  of  curvature  at  the  trailing  edge  (Y^  =  5). 

4.  The  node  spacing  between  nodes  on  the  airfoil  and  those  just  off  it  (10“^). 

5.  The  rate  at  which  node  spacing  grows  close  to  the  airfoil  (1.1). 

6.  The  distance  to  extrude  the  inner  block  (0.1). 

7.  The  half-length  of  a  side  of  the  outer  square  boundary  {R  =  20).  (We  choose 
the  half-length  so  that  it  is  comparable  to  the  radius  of  the  TRANSOM  grid.) 

8.  The  node  spacing  on  the  outer  square  boundary  (2). 

9.  The  boundary  decay  factor  for  the  unstructured  outer  domain.  Values  near  1 
mean  that  the  cell  sizes  change  less  rapidly  (0.99). 

The  grid  is  constructed  as  follows: 

1.  Generate  the  inner  domain  as  if  making  a  TRANSOM  grid,  but  only  extrude 
over  the  size  of  the  inner  grid. 

2.  Make  four  edges  to  form  the  square  outer  boundary  and  populate  them  with 
nodes  using  the  appropriate  spacing. 

3.  Generate  an  unstructured  domain  between  the  outer  boundary  and  the  outer 
edge  of  the  inner  domain  using  the  specified  boundary  decay  factor. 

4.  Use  translational  extrusion  to  convert  the  two-dimensional  grid  into  a  three- 
dimensional  grid  two  cells  deep  (ANSYS  GFX  requires  a  three-dimensional 
grid) . 

5.  Set  boundary  conditions  on  each  exterior  domain. 

6.  Save  the  grid  in  ANSYS  GFX  format. 

An  example  of  a  grid  for  ANSYS  GFX  is  shown  in  Figure  3. 
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Figure  3:  An  example  of  a  grid  generated 
the  left,  a  close-up  near  the  airfoil  on  the 
and  trailing  edges  than  the  grids  actually 


for  ANSYS  CFX.  The  full  grid  is  shown  on 
right.  This  grid  is  coarser  near  the  leading 
used. 


NACA  66  (DTMB  mod):  t  =  0.06  c  =  0 


Figure  4:  The  DTMB  modihed  NACA  66  airfoil  with  zero  cambre 
and  thickness  to  chord  ratio  of  0.06. 

3  Sensitivity  to  grid  parameters 


The  study  of  sensitivity  to  grid  parameters  used  the  DTMB  modified  NACA  66  airfoil 
with  thickness  to  chord  ratio  of  0.06,  roughly  typical  of  a  section  about  halfway  up  a 
propeller  blade,  and  with  no  cambre:  see  Figure  4. 

The  flow  past  the  airfoil  was  calculated  at  an  angle  of  attack  of  4°,  large  enough 
to  generate  a  significant  pressure  peak  with  a  value  of  the  pressure  coefficient,  Cp, 
near  —6.  For  all  calculations  the  Reynolds  number  based  on  the  airfoil  chord  was 
Rf.  =  10’^.  For  the  calculations  with  ANSYS  CFX  the  SST  turbulence  model  [8]  was 
used,  while  for  TRANSOM  the  Spalart-Allmaras  turbulence  model  [9]  was  used. 
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Pressure:  CFX 


Pressure:  TRANSOM 


Figure  5:  Pressure  coefficient  near  the  leading  edge  at  different  values 
ofs;  ANSYS  CFX  on  the  left,  TRANSOM  on  the  right. 


3.1  Sensitivity  to  waii  spacing 


The  separation  of  the  wall  from  the  nodes  closest  to  it  is  nsually  measured  using  the 
parameter  y~^: 


+  _  “1/ 
y  =  — ; 

V 


5);  2  _ 

u  =  v 


dv 

dy 


(3) 


where  v  is  the  kinematic  viscosity,  v  is  the  component  of  the  velocity  parallel  to  the 
wall  and  y  is  the  distance  to  the  wall.  Because  the  value  of  y^  depends  on  the  velocity, 
it  has  a  peak  near  the  leading  edge  where  both  the  velocity  and  pressure  have  peaks. 
Therefore  the  values  of  y^  near  the  pressure  peak  can  be  quite  high  relative  to  the 
values  over  the  main  portion  of  the  airfoil.  Is  it  necessary  to  choose  the  wall-spacing 
based  on  the  high  values  of  y^  in  the  pressure  peak? 


The  wall  spacing  for  the  initial  grid  generated  for  ANSYS  CFX  was  chosen  so  that  the 
largest  value  of  y^  was  about  11.  The  ANSYS  CFX  scalable  wall  functions  effectively 
set  the  minimum  allowed  value  of  t/’*'  to  11,  so  we  expect  to  see  little  change  if  the  wall 
spacing  is  smaller  than  this.  A  series  of  grids  were  then  generated  in  which  the  wall 
spacing  was  smaller  or  larger  than  the  initial  value  by  a  factor,  s.  A  similar  procedure 
was  used  for  the  TRANSOM  grids  except  that  the  wall  spacing  of  the  initial  grid  was 
set  so  that  the  largest  value  of  y^  was  1  (TRANSOM  does  not  use  wall  functions 
so  the  nodes  must  be  closer  to  the  wall  than  for  ANSYS  CFX).  Thus,  if  s  =  2,  the 
largest  value  of  y^  for  the  grid  generated  for  ANSYS  CFX  is  about  22,  while  for  the 
TRANSOM  grid  it  is  about  2.  The  value  s  =  1  represents  a  wall  spacing  that,  a 
priori,  would  be  considered  optimal. 


For  both  solvers  the  spacing  of  the  nodes  at  the  leading  edge  was  set  so  that  Ni^.  =  10. 


Figure  5  shows  the  pressure  coefficient  predicted  by  both  solvers  near  the  leading 
edge  as  a  function  of  the  chordwise  distance  (the  leading  edge  is  0;  the  trailing  edge 
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Skin  Friction:  CFX 


Skin  Friction:  TRANSOM 


0  0.002  0.004  0.006  0.008  0.01  0  0.002  0.004  0.006  0.008  0.01 

Chord  Chord 

Figure  6:  Skin  friction  coefficient  near  the  hading  edge  at  different 
vaiues  of  s;  ANSYS  CFX  on  the  ieft,  TRANSOM  on  the  right. 


Pressure:  CFX  and  TRANSOM 


Figure  7:  Comparison  of  ANSYS  CFX 
and  TRANSOM  predictions  of  the  pres¬ 
sure  coefficient  near  the  hading  edge 
with  s  =  1. 


Skin  Friction:  CFX  and  TRANSOM 


Chord 

Figure  8:  Comparison  of  ANSYS  CFX 
and  TRANSOM  predictions  of  the  skin 
friction  coefficient  near  the  hading  edge 
with  s  =  0.5. 


is  1).  The  ANSYS  CFX  predictions  vary  very  little  for  y~^  as  high  as  350,  while  the 
TRANSOM  predictions  are  jnst  beginning  to  degenerate  when  y+  reaches  32.  Clearly 
the  pressnre  is  insensitive  to  the  wall  spacing.  As  the  pressnre  distribntion  is  almost 
entirely  a  potential  flow  effect,  this  is  what  shonld  be  expected. 

Fignre  6  shows  the  skin  friction  coefficient,  C/,  near  the  leading  edge.  In  this  case 
there  is  a  clear  dependence  on  wall  spacing  for  both  ANSYS  CFX  and  TRANSOM. 
With  ANSYS  CFX,  the  dependence  continnes  even  to  spacings  somewhat  smaller 
than  s  =  1.  The  TRANSOM  resnlts  have  converged  at  s  =  1  bnt  begin  to  diverge 
for  larger  spacings. 

Fignre  7  compares  the  pressnre  and  skin  friction  distribntions  when  s  =  1  predicted 
by  ANSYS  CFX  and  TRANSOM.  There  is  very  good  agreement.  Similarly,  Fignre  8 
compares  the  skin  friction  distribntions  when  s  =  0.5.  The  agreement  is  good  except 
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Figure  9:  Variation  of  lift  (left)  and  drag  (right)  coefficients  with  . 


for  a  small  region  just  aft  of  the  pressure  peak  on  the  suction  side  where  the  ANSYS 
CFX  prediction  is  smaller  than  that  of  TRANSOM.  Therefore,  as  should  be  expected, 
when  each  solver  is  used  with  near-optimal  node  spacing,  they  produce  very  similar 
results. 

It  is  important  to  note  that  when  s  >  1,  the  region  over  which  the  skin  friction  is  not 
well  predicted  is  very  small;  consequently,  the  effect  on  the  lift  and  drag  of  the  airfoil 
is  also  small:  see  Figure  9.  Therefore,  when  gridding  a  propeller,  it  will  be  sufficient 
to  choose  the  near  wall  spacing  so  that  y~^  is  close  to  1  (TRANSOM)  or  11  (ANSYS 
CFX)  over  the  main  portion  of  the  airfoil;  if  it  exceeds  these  values  near  the  pressure 
peak,  there  will  be  little  effect  on  the  predicted  propeller  performance. 


3.2  Sensitivity  to  node  spacing  at  the  ieading  edge 

To  test  the  sensitivity  to  the  spacing  of  the  nodes  near  the  leading  edge,  grids  were 
generated  for  which  Ni^,  the  number  of  nodes  within  one  leading  edge  radius  of 
curvature,  was  3,  5,  10,  15  or  20.  Since  the  radius  of  curvature  at  the  leading  edge  is 
1.63  X  10“^,  these  value  of  Nie  correspond  to  node  spacings  of  5.43  x  10“^,  3.26  x  lO”"^, 
1.63  X  10“^,  1.09  X  lO”"^  and  8.15  x  10“®.  Figure  10  shows  the  nodes  near  the  leading 
edge  (the  same  for  both  ANSYS  CFX  and  TRANSOM  grids)  when  Ni^  is  3. 

In  all  cases  the  wall  spacing  was  chosen  so  that  y~^  has  a  maximum  value  of  1  for 
TRANSOM  and  11  for  ANSYS  CFX. 

Figure  11  shows  the  pressure  coefficient  near  the  leading  edge  as  a  function  of  the 
chordwise  distance.  For  both  ANSYS  CFX  and  TRANSOM  the  pressure  is  well- 
predicted  at  Nie  =  5.  There  is  only  a  small  change  in  the  pressure  peak  when  Ni^  is 
reduced  to  3. 

Figure  12  shows  the  skin  friction  coefficient,  Cf,  near  the  leading  edge.  In  this  case 
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Leading  edge  nodes:  N|g  =  3 


Figure  10:  Nodes  near  the  leading  edge  when  Nie  =  3. 


Pressure  near  the  LE:  CFX 


Pressure  near  the  LE:  TRANSOM 


Chord 


Figure  11:  Pressure  coethcient  near  the  leading  edge  at  different  values  of 
Nie,  the  number  of  nodes  relative  to  the  radius  of  curvature  leading  edge; 
ANSYS  CFX  on  the  left,  TRANSOM  on  the  right. 


Skin  Friction  near  the  LE:  CFX 


Skin  Friction  near  the  LE:  TRANSOM 


Chord 


Chord 


Figure  12:  Skin  friction  coefficient  near  the  leading  edge  at  different  values 
of  Nie,  the  number  of  nodes  relative  to  the  radius  of  curvature  leading  edge; 
ANSYS  CFX  on  the  left,  TRANSOM  on  the  right. 
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Pressure  near  the  LE:  CFX 


Pressure  near  the  LE:  TRANSOM 


Figure  13:  Pressure  coefEcient  near  the  leading  edge  for  different  grid 
sizes;  ANSYS  CFX  on  the  left,  TRANSOM  on  the  right. 


Skin  Friction  near  the  LE:  CFX 


Skin  Friction  near  the  LE:  TRANSOM 


Figure  14:  Skin  friction  coefficient  near  the  leading  edge  for  different 
grid  sizes;  ANSYS  CFX  on  the  left,  TRANSOM  on  the  right. 


the  ANSYS  CFX  predictions  show  a  clear  dependence  on  A/e,  thongh  not  enongh  to 
alter  the  lift  or  drag  signihcantly.  The  TRANSOM  results  are  nearly  independent  of 
Nle. 

3.3  Sensitivity  to  the  size  of  the  grid 

To  test  the  sensitivity  to  the  overall  size  of  the  grid,  R,  grids  were  generated  for  which 
R  was  10,  20,  30  or  40.  In  all  cases  the  wall  spacing  was  chosen  so  that  has  a 
maximum  value  of  1  for  TRANSOM  and  11  for  ANSYS  CFX.  For  both  solvers  the 
spacing  of  the  nodes  at  the  leading  edge  was  set  so  that  A/e  =  10. 

Figure  13  shows  the  pressure  coefficient  near  the  leading  edge  as  a  function  of  the 
chordwise  distance.  For  both  ANSYS  CFX  and  TRANSOM  there  is  little  variation 
with  R,  the  largest  deviation  occurring,  as  expected,  when  R  =  10.  Figure  14  shows 
that  the  skin-friction  coefficient  near  the  leading  edge  is  likewise  insensitive  to  R. 
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Skin  Friction  near  the  LE:  CFX 


Figure  15:  Skin  friction  coefficient  near  the 
leading  edge  for  different  levels  of  turbu¬ 
lence  in  the  free  stream. 

3.4  Sensitivity  of  ANSYS  CFX  to  the  size  of  the 
eiements  near  the  outer  boundary 

To  test  the  sensitivity  to  the  overall  size  of  the  grid,  R,  grids  were  generated  for  which 
the  onter  node  spacing  was  0.5,  1,  2  or  3.  The  effect  on  the  pressnre  distribntion  and 
the  skin-friction  distribntion  was  so  small  as  to  be  indiscernible  when  plotted. 

3.5  Sensitivity  of  ANSYS  CFX  to  the  turbuience 
intensity 

At  inlet  bonndaries,  ANSYS  CFX  reqnires  that  the  valnes  of  tnrbulence  variables  be 
set.  This  sets  the  level  of  tnrbnlence  in  the  free  stream.  In  all  the  rnns  described  here 
this  was  done  by  specifying  the  tnrbnlence  intensity,  I,  and  the  viscosity  ratio  jJLt/ h- 
Three  preset  options  are  provided  called  low  (/  =  1%,  Ht/ h  =  1))  medinm  (J  =  5%, 
git/ b-  =  10)  and  high  (J  =  10%,  git/ gi  =  100).  The  effect  of  these  settings  on  the  skin 
friction  coefficient  near  the  leading  edge  is  shown  in  Fignre  15. 

The  level  of  tnrbnlence  has  a  clear  effect  on  the  skin  friction  jnst  aft  of  the  pressnre 
peak  on  the  snction  side  of  the  blade,  the  region  where  there  was  a  difference  between 
ANSYS  CFX  and  TRANSOM.  Once  again,  the  difference  is  not  snfficiently  large  that 
there  will  be  an  appreciable  effect  on  the  lift  or  drag. 

The  level  of  tnrbnlence  had  no  discernible  effect  on  the  pressnre  distribntion. 
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4  Predictions  of  maximum  pressure 


To  put  the  size  of  the  variations  with  the  grid  parameters  into  perspective  relative 
to  the  variations  one  hnds  between  solvers,  the  maximum  pressure  on  a  series  of 
airfoils  was  calculated  for  different  cambres,  thicknesses  and  angles  of  attack  using 
both  ANSYS  CFX  and  TRANSOM.  The  results  for  the  DTMB  modihed  NACA  66 
airfoil  with  cambre  of  2%  of  chord  are  plotted  as  cavitation  buckets  in  Figure  16.  Each 
curve  shows  the  relationship  between  the  angle  of  attack  and  the  minimum  value  of 
the  pressure  coefficient,  Cp,  on  the  airfoil.  The  different  curves  represent  different 
thickness  to  cambre  ratios. 

Relative  to  TRANSOM,  ANSYS  CFX  tends  to  predict  slightly  higher  peak  pressures 
for  positive  angles  of  attack,  slightly  lower  for  negative  angles  of  attack.  This  results 
in  a  small  upward  shift  of  the  buckets  relative  to  those  predicted  by  TRANSOM.  The 
same  effect  can  be  seen  in  Figure  7.  The  reason  for  this  discrepancy  is  not  understood. 
The  effect  exceeds  the  effects  shown  by  variations  in  the  grid  (for  example,  compare 
Figures  5  and  11  with  Figure  7),  so  it  is  believed  to  be  caused  by  the  inherent 
differences  in  the  solvers. 

When  the  thickness  t  =  0.02  the  TRANSOM  curve  shows  the  minimum  pressure 
dropping  precipitously  once  the  angle  of  attack  gets  too  large.  This  is  due  to  separa¬ 
tion  beginning  near  the  leading  edge  causing  the  airfoil  to  stall.  The  lack  of  similar 
behaviour  in  the  ANSYS  CFX  shows  that  the  airfoil  has  not  yet  stalled  even  at  sig- 
nihcantly  larger  angles  of  attack.  This  discrepancy  is  almost  certainly  due  to  the 
different  turbulence  models  used:  SST  for  ANSYS  CFX  and  Spalart-Allmaras  for 
TRANSOM. 

These  results  suggest  that  errors  associated  with  grids  constructed  using  the  nominal 
parameters  given  in  Sections  2.1  and  2.2  will  be  smaller  than  the  modelling  errors 
associated  with  the  different  solvers  and  turbulence  models. 

5  Concluding  remarks 


The  grid  sensitivity  of  the  RANS  solvers  ANSYS  CFX  and  TRANSOM  has  been 
assessed  for  two-dimensional  flows  around  propeller  blade  sections  at  small  angles  of 
attack.  Both  solvers  show  little  sensitivity  to  the  grid  parameters  tested.  These  results 
suggest  that  accurate  predictions  of  the  leading  edge  pressure  peaks  on  propeller 
blades  can  be  attained  with  ANSYS  CFX  provided  that  the  values  of  the  near 
wall  nodes  is  within  the  log-layer  of  the  boundary  layer  {y~^  ~  300)  and  the  leading 
edge  node  density  is  such  that  at  least  hve  nodes  are  contained  within  one  radius  of 
curvature.  The  distance  to  the  outer  boundary  of  the  grid  should  exceed  10  chord 
lengths.  Similarly,  TRANSOM  will  require  y~^  values  not  exceeding  about  16.  It  is 
not  so  clear  whether  these  restrictions  will  be  adequate  to  predict  the  complicated 
vortex  production  near  the  propeller  tip. 
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t  =  0.20  t  =  0.18  t  =  0.16  t  =  0.14  t  =  0.12 
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Figure  16:  Cavitation  buckets  for  the  DTMB  modified  NACA  66  airfoil 
with  2%  cambre  ratio  as  calculated  by  TRANSOM  and  ANSYS  CFX. 
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List  of  symbols 


/i  Viscosity. 

fit  Eddy  viscosity. 

u  Kinematic  viscosity. 

Cf  The  skin  friction  coefficient. 

Cp  The  pressure  coefficient. 

I  Turbulence  intensity. 

Nie  The  number  of  nodes  per  radius  of  curvature  at  the  leading  edge. 

Nte  The  number  of  nodes  per  radius  of  curvature  at  the  leading  edge. 

R  The  outer  radius  of  a  TRANSOM  grid.  The  half  length  of  the  outer  square  of 
a  grid  used  by  ANSYS  CFX. 

s  A  factor  used  to  indicate  the  relative  spacing  of  the  near-wall  nodes  in 
different  grids. 

sie  Spacing  of  the  nodes  on  the  airfoil  surface  at  the  leading  edge. 

Sm  Mean  spacing  of  the  nodes  on  the  airfoil  surface. 

Sw  The  spacing  between  nodes  on  the  airfoil  and  those  just  off  it. 

t  Thickness  to  chord  ratio  of  an  airfoil. 

u*  Shear  velocity:  see  Equation  (3). 

w  The  growth  rate  of  the  extrusion  of  the  structured  block. 

Non-dimensional  distance  to  the  airfoil  surface:  see  Equation  (3). 
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